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Abstract. In this paper we are concerned with the monodromy of Picard- 
Fuchs differential equations associated with one-parameter families of Calabi- 
Yau threefolds. Our results show that in the hypergeometric cases the matrix 
representations of monodromy relative to the Frobenius bases can be expressed 
in terms of the geometric invariants of the underlying Calabi-Yau threefolds. 
This phenomenon is also verified numerically for other families of Calabi-Yau 
threefolds in the paper. Furthermore, we discover that under a suitable change 
of bases the monodromy groups are contained in certain congruence subgroups 
of Sp(4, Z) of finite index and whose levels arc related to the geometric invari- 
ants of the Calabi-Yau threefolds. 



1. Introduction 

Let M z be a family of Calabi-Yau n-folds parameterized by a complex variable 
z G P 1 (C), and lo z be the unique holomorphic differential n-form on M z (up to 
a scalar). Then the standard theory of Gauss-Manin connections asserts that the 
periods 



L 



L0 Z 

'7* 

satisfy certain linear differential equations, called the Picard-Fuchs differential equa- 
tions, where "f z are r-cycles on M z . 

When n= 1, Calabi-Yau onefolds are just elliptic curves. A classical example of 
Picard-Fuchs differential equations is 

(1) (1 - z)9 2 f - z9f - Z -f = 0, 9 = zd/dz, 

satisfied by the periods 

dx 

■\/x(x — l)(x — z) 

of the family of elliptic curves E z : y 2 = x(x — l)(x — z). 

When n = 2, Calabi-Yau manifolds are either 2-dimcnsional complex tori or K3 
surfaces. When the Picard number of a one-parameter family of K3 surfaces is 1 9, 
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the Picard-Fuchs differential equation has order 3. One of the simplest examples is 

x\ + x\ + x\ + x\ — Z~ 1 XiX 2 X3X4 = C P 3 , 

whose Picard-Fuchs differential operator is 

(2) 3 - 42(46* +l)(40 + 2)(40 + 3). 
Another well-known example is 

(3) (1 - 34z + z 2 )6 3 + (3z 2 - blz)6 2 + (3z 2 - 27z)9 + {z 1 - 5z), 
which is the Picard-Fuchs differential operator for the family of K3 surfaces 

1 - (1 - XY)Z - zXYZ{\ - X)(l -Y)(l-Z) = 0. 

(See [7].) This differential equation appeared in Apery's proof of irrationality of 
C(3). (See 0.) 

When n = 3 and Calabi-Yau threefolds have the Hodge number h 2 ' 1 equal to 1, 
the Picard-Fuchs differential equations have order 4. One of the most well-known 
examples of such Calabi-Yau threefolds is the quintic threefold 

x\ + x\ + x\ + x\ + x\ — z^ 1 11X2X3X^X5 = C P 4 . 

In [9], it is shown that the Picard-Fuchs differential operator for this family of 
Calabi-Yau threefolds is 

(4) 4 - 5z(50 + 1)(50 + 2) (50 + 3) (50 + 4). 

(Actually, it is the mirror partner of the quintic Calabi-Yau threefolds that has 
Hodge number ft 2 ' 1 = 1 and hence the Picard-Fuch differential equation is of or- 
der 4. But the mirror pair of Calabi-Yau threefolds share the same "principle 
periods" . This means that the Picard-Fuchs differential equation of the original 
quintic Calabi-Yau threefold of order 204 contains the above order 4 equation as a 
factor and the factors corresponding to the remaining 200 "semiperiods" . 

In this article we are concerned with the monodromy aspect of the Picard-Fuchs 
differential equations. Let 

L : r n (z)9 n + r n _ 1 ^)0"- 1 + • ■ ■ + r Q (z), r t e C(z), 

be a differential operator with regular singularities. Let zq be a singular point 
and S be the solution space of L at Zq. Then analytic continuation along a closed 
curve 7 circling zo gives rise to an automorphism of S, called monodromy. If a basis 
{/1, ...,/„} of S is chosen, then we have a matrix representation of the monodromy. 
Suppose that fa becomes a,i/i + ■ • • + ai„f n after completing the loop 7, that is, if 

(an ■■■ ai n \ ( fi 



\fnj \^nl • ■ ■ Q>nn / \fn; 

then the matrix representation of the monodromy along 7 relative to the basis 
{fi} is the matrix (a^). The group of all such matrices are referred to as the 
monodromy group relative to the basis {fi} of the differential equation. Clearly, 
two different choices of bases may result in two different matrix representations 
for the same monodromy. However, it is easily seen that they are connected by 
conjugation by the matrix of basis change. Thus, the monodromy group is denned 
up to conjugation. In the subsequent discussions, for the ease of exposition, we 
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may often drop the phrase "up to conjugation" about the monodromy groups, 
when there is no danger of ambiguities. 

It is known that for one-parameter families of Calabi-Yau varieties of dimension 
one and two (i.e., elliptic curves and K3 surfaces, respectively), the monodromy 
groups are very often congruence subgroups of SX(2, R). For instance, the mon- 
odromy group of (P) is T(2), while those of ((2]) and ((3]) are r (2)-|-u;2 and T (6)+ujq, 
respectively, where u>d denotes the Atkin-Lehner involution. (Technically speaking, 
the monodromy groups of @ and are subgroups of SX(3,R) since the order 
of the differential equations is 3. But because © and §5§ are symmetric squares 
of second-order differential equations, we may describe the monodromy in terms 
of the second-order ones.) Moreover, suppose that yo(z) = 1 + • • • is the unique 
holomorphic solution at z — and y\{z) = yo(z) logz + g{z) is the solution with 
logarithmic singularity. Set r = cyi(z) /yo(z) for a suitable complex number c. 
Then z, as a function of r, becomes a modular function, and yo(z(T)) becomes a 
modular form of weight 1 for the order 2 cases and of weight 2 for the order 3 cases. 
For example, a classical result going back to Jacobi states that 



or equivalently, that the modular form y(r) = 0|, as a function of z{t) — B\jB\, 
satisfies ((TJ). Here 2 F\ denotes the Gauss hypergeometric function. 

In this paper we will address the monodromy problem for Calabi-Yau threefolds. 
At first, given the experience with the elliptic curve and K3 surface cases, one may 
be tempted to guess that the monodromy group of such a differential equation will 
be the symmetric cube of some congruence subgroup of SX(2, R). After all, there 
is a result by Stiller [23 (see also [27 ) asserting that if t(r) is a non-constant mod- 
ular function and F(t) is a modular form of weight k on a subgroup of SL(2,R) 
commensurable with <SX(2, Z), then F, tF, . . . , r k F, as functions of t, are solutions 
of a (k + l)-st order linear differential equation with algebraic functions of t as 
coefficients. However, this is not the case in general. A quick way to see this is 
that the coefficients of the symmetric cube of a second order differential equation 
y" + rx(t)y' + r (t)y = is completely determined by r\ and r , but the coeffi- 
cients of the Picard-Fuchs differential equations, including (j4|), do not satisfy the 
required relations. (The exact relations can be computed using Maple's command 
symmetric_power.) Nevertheless, in the subsequent discussion we will show that, 
with a suitable choice of bases, the monodromy groups for Calabi-Yau threefolds 
are contained in certain congruence subgroups of Sp(4, Z) whose levels are somehow 
described in terms of the geometric invariants of the manifolds in question. This 
is proved rigorously for the hypergeometric cases and verified numerically for other 
(e.g., non- hypergeometric) cases. Furthermore, our computation in the hypergeo- 
metric cases shows that the matrix representation of the monodromy around the 
finite singular point (different from the origin) relative to the Frobenius basis at 
the origin can be expressed completely using the geometric invariants of the asso- 
ciated Calabi-Yau threefolds. This phenomenon is also verified numerically in the 
non-hypergeometric cases. Although it is highly expected that geometric invariants 




where 
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will enter into the picture, in reality, geometry will dominate the entire picture in 
the sense that every entry of the matrix is expressed exclusively in terms of the 
geometric invariants. 

The monodromy problem in general has been addressed by a number of authors. 
Papers relevant to our consideration include [6], [9], [11], [16] , and [26] ■ to name 
a few. In [6] , Beukers and Heckman studied monodromy groups for the hypergeo- 
metric functions n F n —i. They showed that the Zariski closure of the monodromy 
groups of (|3]) is Sp(4, C). The same is true for other Picard-Fuchs differential equa- 
tions for Calabi-Yau threefolds that are hypergeometric. In [5], Candelas et al. 
obtained precise matrix representations of monodromy for ((4]). Then Klemm and 
Theisen [16] applied the same method as that of Candelas et al. to deduce mon- 
odromy groups for three other hypergeometric cases. In |llj Doran and Morgan 
determined the monodromy groups for all the hypergeometric cases. Their matrix 
representations also involve geometric invariants of the Calabi-Yau threefolds. For 
Picard-Fuchs differential equations of non-hypergeometric type, there is not much 
known in literature. In |26] van Enckevort and van Straten computed the mon- 
odromy matrices numerically for a large class of differential equations. In many 
cases, they are able to find bases such that the monodromy matrices have rational 
entries. We will discuss the above results in more detail in Sections 3-5. 

Our motivations of this paper may be formulated as follows. Modular functions 
and modular forms have been extensively investigated over the years, and there are 
great body of literatures on these subjects. As we illustrated above, the monodromy 
groups of Picard-Fuchs differential equations for families of elliptic curves and K2> 
surfaces are congruence subgroups of SL(2,M). This modularity property can be 
used to study properties of the differential equations and the associated manifolds. 
For instance, in [18] Lian and Yau gave a uniform proof of the integrality of Fourier 
coefficients of the mirror maps for several families of K3 surfaces using the fact 
that the monodromy groups are congruence subgroups of SL2{M). For such an 
application, it is important to express monodromy groups in a proper way so that 
properties of the associated differential equations can be more easily discussed and 
obtained. Thus, the main motivation of our investigation is to find a good rep- 
resentation for monodromy groups from which further properties of Picard-Fuchs 
differential equations for Calabi-Yau threefolds can be derived. 

The terminology "modularity" has been used for many different things. One 
aspect of the modularity that we would like to address is the modularity question 
of the Galois representations attached to Calabi-Yau threefolds, assuming that 
Calabi-Yau threefolds in question are defined over Q. Let X be a Calabi-Yau 
threefold defined over Q. We consider the L-series associated to the third etale 
cohomology group of X. It is expected that the L-series should be determined 
by some modular (automorphic) forms. The examples of Calabi-Yau threefolds we 
treat in this paper are those with the third Betti number equal to 4. It appears that 
Calabi-Yau threefolds with this property are rather scarce. Batyrev and Straten 
[4] considered 13 examples of Calabi-Yau threefolds with Picard number h 1 ' 1 — 
1. Then their mirror partners will fulfill this requirement. (We note that more 
examples of such Calabi-Yau threefolds were found by Borcea [8].) All these 13 
Calabi-Yau threefolds are defined as complete intersections of hypersurfaces in 
weighted projective spaces, and they have defining equations defined over Q. 
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To address the modularity, we ought to have some "modular groups" , and this 
paper offers candidates for appropriate modular groups via the monodromy group 
of the associated Picard-Fuchs differential equation (of order 4). In these cases, we 
expect that modular forms of more variables, e.g., Siegel modular forms associated 
to the modular groups for our congruence subgroups would enter the scene. 

In general, the third Betti numbers of Calabi-Yau threefolds are rather large, 
and consequently, the dimension of the associated Galois representations would be 
rather high. To remedy this situation, we first decompose Calabi-Yau threefolds 
into motives, and then consider the motivic Galois representations and their mod- 
ularity. Especially, when the principal motives (e.g., the motives that are invariant 
under the mirror maps) are of dimension 4, the modularity question for such mo- 
tives should be accessible using the method developed for the examples discussed 
in this paper. 

The modularity questions will be treated in subsequent papers. 

2. Statements of results 

To state our first result, let us recall that among all the Picard-Fuchs differential 
equations for Calabi-Yau threefolds, there are 14 equations that are hypergeometric 
of the form 

6» 4 - Cz(9 + A) (6 + 1 - A)(6 + B)(9 + 1 - B). 
Their geometric descriptions and references are given in the following Table 1. 



# 


A 


B 


C 


Description 


H 3 


c 2 -H 


C3 


Ref 


1 


1/5 


2/5 


3125 


X(5) C P 4 


5 


50 


-200 


m 


2 


1/10 


3/10 


8 • 10 5 


X(10)cP 4 (l, 1,1,2,5) 


1 


34 


-288 


\M 


3 


1/2 


1/2 


256 


X(2,2,2,2) C P 7 


16 


64 


-128 


M 


4 


1/3 


1/3 


729 


X{3,3) C P 5 


9 


54 


-144 


M 


5 


1/3 


1/2 


432 


X(2,2,3) C P 6 


12 


60 


-144 


m 


6 


1/4 


1/2 


1024 


X(2,4) c P 5 


8 


56 


-176 


m 


7 


1/8 


3/8 


65536 


X(8) C P 4 (l,l,l,l,4) 


2 


44 


-296 


[20] 


8 


1/6 


1/3 


11664 


X(6)cP 4 (l,l,l,l,2) 


3 


42 


-204 


[20] 


9 


1/12 


5/12 


12 6 


X{2, 12) CP 5 (1,1,1, 1,4,6) 


1 


46 


-484 


[IT] 


10 


1/4 


1/4 


4096 


JC(4,4) C P 5 (l,l,l,l,2,2) 


4 


40 


-144 


PI] 


11 


1/4 


1/3 


1728 


X(4,6) C P 5 (l,l,l,2,2,3) 


6 


48 


-156 


PI] 


12 


1/6 


1/4 


27648 


X(3,4) cP 5 (l,l,l,l,l,2) 


2 


32 


-156 


PI] 


13 


1/6 


1/6 


2 8 -3 6 


X{6,6) c P 5 (l,l,2,2,3,3) 


1 


22 


-120 


PU 


14 


1/6 


1/2 


6912 


Y(2,6) cP 5 (l, 1,1, 1,1,3) 


4 


52 


-256 


PI] 



Some comments might be in order for the notations in the table. We employ the no- 
tations of van Enckevort and van Straten [26]. X(d\, d%, ■ ■ ■ , dk) C P n (wo, . . . , w n ) 
stands for a complete intersection of k hypersurfaces of degree d\ , . . . , dk in the 
weighted projective space with weight (u> , • • • , w n ). For instance, X(3, 3) C P 5 is 
a complete intersection of two cubics in the ordinary projective 5-space P 5 defined 

by 

{ if + Y 2 3 + Y 3 3 - 30Y 4 Y 5 r 6 - } n { -30YiY 2 Y 3 + Y 4 3 + Y 5 3 + Y 6 3 = }. 
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Slightly more generally, Y(4, 4) C P 5 (l, 1, 2, 1, 1, 2) denotes a complete intersection 
of two quartics in the weighted projective 5-space P 5 (l, 1, 2, 1, 1, 2) and may be 
defined by the equations 

{ y 4 + y 4 + y 2 - 44>y a y b y 9 = o } n { -^y x y 2 y z + y 4 4 + y 5 4 + y 6 2 = o }. 

We note that all these examples of Calabi-Yau threefolds M have the Picard number 
h l,1 {M) = 1. Let 0(H) be the ample generator of the Picard group Pic(M) ~ Z. 
The basic invariants for such a Calabi-Yau threefold M are the degree d := H 3 , 
the second Chern number ci ■ H and the Euler number C3 (the Euler characteristic 
of Al). The equations are numbered in the same way as in pQ. 

In [9], using analytic properties of hypergeometric functions, Candelas et al. 
proved that with respect to a certain basis, the monodromy matrices around z = 
and z = 1/3125 for the quintic threefold case (Equation 1 from Table 1) are 



/51 


90 


-25 


°\ 





1 








100 


175 


-49 





V-75 


-125 


35 


V 



and 



1 






1 



V 



Applying the 
also obtained 



respectively. (Note that these two matrices are both in Sp(4, Z).) 
same method as that of Candelas et al., Klemm and Theisen [16] 
the monodromy of the one-parameter families of Calabi-Yau threefolds for Equa- 
tions 2, 7, and 8. Presumably, their method should also work for several other 
hypergeometric cases. However, the method fails when the indicial equation of the 
singularity 00 has repeated roots. To be more precise, it does not work for Equa- 
tions 3-6, 10, 13 and 14. Moreover, the method uses the explicit knowledge that 
the singular point z = 1/C is of conifold type. (Note that in geometric terms, a 
conical singularity is a regular singular point whose neighborhood looks like a cone 
with a certain base. For instance, a 3-dimcnsional conifold singularity is locally 
isomorphic to XY - ZT = or equivalently, to X 2 + Y 2 + Z 2 + T 2 = 0. Reflecting 
to the Picard- Fuchs differential equations, this means that the local monodromy is 
unipotent of index 1.) Thus, it can not be applied immediately to study monodromy 
of general hypergeometric differential equations. 

In [11] Doran and Morgan proved that if the characteristic polynomial of the 
monodromy around 00 is 

x 4 + (k - A)x 3 + (6 - 2k + d)x 2 + (k - A)x + 1, 

then there is a basis such that the monodromy matrices around z 
are 



and z = 1/C 



(5) 



(\ 








1 



/ 1 



and 



V 






V 



respectively. It turns out that these numbers d and k both have geometric inter- 
pretation. Namely, the number d = H 3 is the degree of the associated threefolds 
and k = ci-HjYl + H 3 /6 is the dimension of the linear system \H\. Doran and 
Morgan's representation has the advantage that the geometric invariants can be 
read off from the matrices directly (although there is no way to extract the Euler 
number C3 from the matrices), but has the disadvantage that the matrices are no 
longer in the symplectic group (in the strict sense). 
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Before we state our Theorem 1, let us recall the definition of Frobenius basis. 
Since the only solution of the indicial equation at z = for each of the cases is 
with multiplicity 4, the monodromy around z = is maximally unipotent. (See 
[20] for more detail.) Then the standard method of Frobenius implies that at z = 
there are four solutions yj, j = 0, . . . ,, with the property that 

• , yi = y logz + g 1 , 



2/0 = 1 + 

(6) 1,2^, 

2/2 = log z + g 1 \ogz 



1 , 3 1 , 2 

92, 2/3=g2/olog z+-gilog z 



where gi are all functions holomorphic and vanishing at z 
these solutions satisfy the relation 



0. 



+ 92 logz + g 3 , 
We remark that 



Vo 


2/3 




2/1 


2/2 


y'o 


2/3 




2/i 


2/2 



and therefore the monodromy matrices relative to the ordered basis {yo, y2, 2/3, 2/1} 
are in Sp(4, C), as predicted by [6]. Now we can present our first theorem. 

Theorem 1. Let 

L:9 4 - Cz(9 + A){6 + 1 - A){6 + B){0 + 1-B) 

be one of the 14 hypergeometric equations, and H 3 , C2-H, and C3 be geometric 
invariants of the associated Calabi-Yau threefolds given in the table above. Let 
yj, j = 0, ... ,3, be the Frobenius basis specified by Then with respect to the 
ordered basis {y^/(2Tri) 3 , y2/ {2iri) 2 , yi/(2i:i), yo}, the monodromy matrices around 
z — and z — l/C are 



(7) 



respectively, where 



(I 


1 


1/2 


l/6\ 





1 


1 


1/2 








1 


1 


v° 








1 J 



and 



fl + a 

-b 


V ~d 



ab/d 
-b 2 /d 
1 

-b 



a 2 jd \ 
~ab/d 


1 - a J 



C3 



rC(3), 



6 = c 2 -ff/24, 



(27Ti) 3 

Remark 1. We remark that by conjugating by the matrix 



H 







6 
d/2 
1 




d/6 + b 



J 



we do recover Doran and Morgan's representation ([5p. Thus, our Theorem 1 stren- 
thens the results of Doran and Morgan [11]. Although the referee suggested that 
Theorem 1 might be a reformulation of the results of Doran and Morgan. However, 
we do not believe that is the case. For one thing, the argument of Doran-Morgan is 
purely based on Linear Algebra. It might be possible to derive our Theorem 1 com- 
bining the results of Doran-Morgan and those of Kontsevich; we will not address 
this question here, but left to future investigations. 

The appearance of the geometric invariants C2, C3, FL and d is not so surprising. 
In [9], it was shown that the conifold period, defined up to a constant as the 



holomorphic solution f(z) — a\{z ~ l/C) + a^{z — l/C) + 



at z 



l/C that 
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appears in the unique solution f(z) log(z 
at z = 1/C, is asymptotically 

H 3 

(8) 



1/C) +g(z) with logarithmic singularity 



3 c 2 - H , 

lo S z + TS — 7 1°S 2 

48m 



C3 



rC(3) 



6(27ri) 3 b 48™ b (2tu) 3 

near z = 0. (See also [E].) Therefore, it is expected that the entries of the mon- 
odromy matrices should contain the invariants. However, it is still quite remarkable 
that the matrix is determined completely by the invariants alone. We have numer- 
ically verified the phenomenon for other families of Calabi-Yau threefolds, and also 
for general differential equations of Calabi-Yau type. (See pQ for the definition of a 
differential equation of Calabi-Yau type. See also Section 5 below.) It appears that 
if the differential equation has at least one singularity with exponents 0, 1, 1, 2, then 
there is always a singularity whose monodromy relative to the Frobenius basis is 
of the form stated in the theorem. Thus, this gives a numerical method to identify 
the possible geometric origin of a differential equation of Calabi-Yau type. 

We emphasize that our proof of Theorem 1 is merely a verification. That is, 
we can prove it, but unfortunately it does not give any geometric insight why the 
matrices are in this special form. 

Acutally, the referee has pointed out that such a geometric interpretation seems 
to exist by Kontsevich. In the framework of "homological mirror symmetry" of 
Kontsevich, the first matrix in Theorem 1 would be the matrix associated to ten- 
soring by the hyperplane line bundle in the bounded derived category of sheaves on 
the Calabi-Yau variety. In general, the matrices in Theorem 1 describe the coho- 
mology action of certain Fourier-Mukai functors. In particular, this explains why 
the matrices are determined by topological invariants of the underlying Calabi-Yau 
manifolds. The paper of van Enckevort and van Straten [26j addressed monodromy 
calculations of fourth order equations of Calabi-Yau type based on homological 
mirror symmetry. The reader is referred to the article |26j for full details about 
geometric interpretations of matrices. We wonder, though, if the Kontsevich's 
results fully explain why there are no "non-geometric" numbers in the second ma- 
trices. To be more precise, here is our question. Since the second matrix M is 
unipotent of rank 1, we know that the rows of M — Id are all scalar multiples of a 
fixed row vector. We probably can deduce from Kontsevich's result that the fourth 
row is (— d, 0, —6, —a), but why the first three rows of M — Id are —a/d, —b/a, and 
times this vector (but not other "non-geometric" scalars? 

Now conjugating the matrices by 



(9) 



/ 



\-d 






d 





1 



d/2 

-b 



1 

-b 
-a) 



we can bring the matrices into the symplectic group Sp(4,Z). The results are 

A 
o 

d 



1 









f 1 








°\ 


1 








and 





1 





1 


d 


1 











1 





-k 


-1 


l ) 












V 



for z = and z = 1/C, respectively, where k = 2b + d/6. Since the monodromy 
group is generated by these two matrices, we see that the group is contained in 
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the congruence subgroup T(d, gcd(d, k)), where the notation T(di,d,2) with d^\d\ 
represents 



r(di,d2) = ^7eSp(4,Z): 7 i 



n 



A 


* 


* 











* 


* 


* 




mod d 





1 












\o 


* 


* 


*) 












(I 


* 


* *\ 











i 


* * 


Z) 


7 










1 













* ij 



mod c?2 



We remark that the entries of the matrices in r(di, ^2) satisfy certain congruence 
relations inferred from the symplecticity of the matrices. To be more explicit, let 
us recall that the symplectic group is characterized by the property that 



7 



^ GSp(2n,C), 



where A, B, C, and D are n x n blocks, if and only if 



7 



D l -B t 
-C* A t 



Thus, for 

/ an ai2 013 014^ 

021 022 «23 a 24 
031 032 O33 a 34 

\a 4 i 042 043 a 44y / 
to be in r(di,d2), the integers a,j should satisfy the implicit conditions 

a 22 a 44 — a 24 a 42 = 1, a 23 = a 14 a 22 — a 12 a 24, a 43 = a 14 a 42 — a 12 a 44 mod dl, 

and 

ai2 = — 043 mod d2- 
We now summarize our finding in the following theorem. 

Theorem 2. Let 

9 4 - Cz(9 + A){6 + 1 - A)(6 + B){6 + 1 - B) 

0, . . . , 3 be the Frobenius 



be one of the 14 hypergeometric equations. Let yj, j 
basis. Then relative to the ordered basis 



2m 



J/o, 



H 3 H 3 



c 2 -H 



-2/0, 



H 3 



2/3 - 



c 2 -H 



24 6(27ri) 3 
i/ie monodromy matrices around z — and z = 1/C are 



(10) 



487ri 



72/1' 



C3 



(2ttz) 



:C(3)l/0, 



(I 


1 





°\ 




n 








°\ 





1 








and 





1 





1 


d 


d 


1 











1 







-jfe 


-1 














V 
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with d — H 3 , k — H 3 /6 + C2-H/12, respectively. They are contained in the congru- 
ence subgroups T(di,d 2 ) for the 14 cases in the table below. 



# 


A B 


d\ d 2 


# 


A B 


di d 2 


1 


1/5 2/5 


5 5 


8 


1/6 1/3 


3 1 


2 


1/10 3/10 


1 1 


9 


1/12 5/12 


1 1 


3 


1/2 1/2 


16 8 


10 


1/4 1/4 


4 4 


4 


1/3 1/3 


9 3 


11 


1/4 1/3 


6 1 


5 


1/3 1/2 


12 1 


12 


1/6 1/4 


2 1 


6 


1/4 1/2 


8 2 


13 


1/6 1/6 


1 1 


7 


1/8 3/8 


2 2 


14 


1/6 1/2 


4 1 



Remark. We remark that what we show in Theorem 2 is merely the fact that the 
monodromy groups are contained in the congruence subgroups T(di, d 2 ). Although 
the congruence subgroups T(di, d 2 ) are of finite index in Sp(4, Z) (sec the appendix 
by Cord Erdenberger for the index formula), the monodromy groups themselves 
may not be so. 

At this moment, we cannot say anything definite about these groups, e.g., their 
finite indexness. In fact, there are two opposing speculations (one by the authors, 
and the other by Zudilin) about these groups. We believe, based on a result (The- 
orem 13.3) of Sullivan [24], that it might be justified to claim that the monodromy 
group is an arithmetic subgroup of the congruence subgroup T(di,d 2 ), and hence 
is of finite index. (Andrey Todorov pointed out to us Sullivan's theorem, though 
we must confess that we do not fully understand the paper of Sullivan.) 

As opposed to our belief, Zudilin has indicated to us via e-mail that a heuristic 
argument suggests that the monodromy groups are too "thin" to be of finite index. 



It would not be of much significance if the hypergeometric equations are the 
only cases where the monodromy groups are contained in congruence subgroups. 
Our numerical computation suggests that there are a number of further examples 
where the monodromy groups continue to be contained in congruence subgroups of 
Sp(4, Z). However, the general picture is not as simple as that for the hypergeo- 
metric cases. 

As mentioned earlier, our numerical data suggest that the Picard-Fuchs differen- 
tial equations for Calabi-Yau threefolds known in literature all have bases relative 
to which the monodromy matrices around the origin and some singular points of 
conifolds take the form described in Theorem 1. Thus, with respect to the 
basis given in Theorem 2, the matrices around the origin and the conifold points 
again have the form (| 10[) . However, with this basis change, the monodromy matrices 
around other singularities may not be in Sp(4, Z), but in Sp(4, Q) instead, although 
the entries still satisfy certain congruence relations. Furthermore, in most cases, 
we can still realize the monodromy groups in congruence subgroups of Sp(4, Z), by 
a suitable conjugation. 
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Example 1. Consider the differential equation 

250 4 - 15z(510 4 + 840 3 + 720 2 + 300 + 5) 

+ 6z 2 (5310 4 + 8286> 3 + 5416* 2 + 1550 + 15) 

- 54z 3 (4230 4 + 2 1 606> 3 + 43990 2 + 37950 + 1170) 

+ 243z 4 (2790 4 + 13680 3 + 227O0 2 + 15860 + 402) - 59O49z 5 (0 + l) 4 . 

In [4] it is shown that this is the Picard-Fuchs differential equation for the Calabi- 
Yau threefolds defined as the complete intersection of three hypersurfaces of degree 
(1, 1, 1) in P 2 x P 2 x P 2 . The invariants are H 3 = 90, c 2 -H = 108, and c 3 = -90. 
There are 6 singularities 0, 1/27, ±i/V27, 5/9, and oo for the differential equation. 
Among them, the local exponents at z = 5/9 are 0,1,3,4 and we find that the 
monodromy around z = 5/9 is the identity. For others, our numerical computation 
shows that relative to the basis in Theorem 2 the monodromy matrices are 

/l 0\ 

10 1 

10' 

\0 1/ 

/-ll 3 1/3 -l\ 

36 -8 -1 3 

-432 108 13 -36 ' 

\ 108 -27 -3 10 / 

From these, we see that the monodromy group is contained in the following group 

{(oy) G Sp(4,Q) : aij e Z V(t,j) ± (1,3), a 13 £ ^Z, 

a-21, a>3i, a 4i , a 32 , a 34 = mod 18, 011,033 = 1 mod 6, 
042 = 0, 022,044 = 1 mod 3| 
Conjugating by 

/3 0\ 

10 

10' 
\0 1/ 

we find that the monodromy group can be brought into the congruence subgroup 
r(6,3). 

Example 2. Consider the differential equation. 

90 4 - 3z(1730 4 + 34O0 3 + 2720 2 + 1020 + 15) 

- 2z 2 (11290 4 + 5O320 3 + 75970 2 + 47730 + 1083) 
+ 2z 3 (8430 4 + 26280 3 + 23530 2 + 6750 + 6) 

- z 4 (2950 4 + 6O80 3 + 4780 2 + 1740 + 26) + z 5 (0 + l) 4 

This is the Picard-Fuchs differential equation for the complete intersection of 7 
hyperplanes with the Grassmannian G(2, 7) with the invariants H 3 = 42, czH = 84, 
and c 3 = —98. (See [3].) The singularities are 0, 3, 00, and the three roots 
zi = 0.01621 . . ., z 2 = -0.2139 . . ., and z 3 = 289.197 ... of z 3 - 289z 2 - 57z + 



i/VTf : 







(I 


1 





°\ 




To = 





1 










90 


90 


1 









y0 


-24 


-1 


V 


/ 


-17 


3 


1/3 


M 






-54 


10 


1 


3 






-972 


162 


19 


54 




V 


162 


-27 


-3 


-8) 





T 



Tl/27 — 



i/V27 
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1. The monodromy around z = 3 is the identity. The others have the matrix 
representations 



T n = 



T — 

± z 2 — 



/I 



42 



1 
1 

42 
-14 








-1 1/ 

-2\ 



/-13 7 1 

28 -13 -2 4 

-196 98 15 -28 

\ 98 -49 -7 15 J 



(l 0\ 

10 1 

10 

\0 1/ 

/ 1 \ 

42 10 9 

-196 1 -42 

\ 1 j 



Thus, the monodromy group is contained in the subgroup T(14, 7). 



3. A GENERAL APPROACH 

Let 

+ Twit/"" 1 ) + • • • + r 1? / + r y - 0, r, e C(z), 

be a linear differential equation with regular singularities. Then the monodromy 
around a singular point zq with respect to the local Frobcnius basis at zq is actually 
very easy to describe, as we shall see in the following discussion. 

Consider the simplest cases where the indicial equation at zq has n distinct roots 
Ai, . . . , A„ such that A^ — Xj g Z for all i ^ j. In this case, the Frobenius basis 
consists of 

Vj(z) = {z~ z o) Xj .fj{z), j = l,...,n, 

where fj(z) are holomorphic near z and have non- vanishing constant terms. It is 
easy to see that the matrix of the monodromy around z with respect to {yj} is 
simply 

/ e 27riAi o ... o \ 

e 2mX * ■■■ 

\ ••• e 27liX " J 

Now assume that the indicial equation at zq has Ai, . . . , with multiplicities 
ei, . . . , e/c, as solutions, where e\ + ■ ■ ■ + e,k = n and Ai — Xj £ Z for all i ^ j. Then 
for each Xj , there are ej linearly independent solutions 

Vjfi = {z- z ) X \f jfi , 

Vj.l = Vjfi log(z - Z ) + {Z~ Z ) X \fj A , 

Vj.2 = ^Vj,0 log 2 (2 - Zq) + {Z- Z ) Xj fj,l \og(z - Zj) + {z - Z ) X] f ja , 



1 

Vj^-l = {z - 2 ) Aj Jjfj;e 3 -l-h\og h (z - Z ), 

where fj^h are holomorphic near z = zq and satisfy fj,o{zo) = 1 and fj,h{zo) = 
for h > 0. Since jj.h are all holomorphic near zq, the analytic continuation along 
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a small closed curve circling zq does not change fj.h- For other factors, circling zo 
once in the counterclockwise direction results in 



(z - z ) 



and 



log(z - z ) 1 — > log(2 - zo) + 2?ri. 
Thus, the behaviors of yj t h under the monodromy around zq are governed by 





U>4 







(27Tt)' 
( e 3" 



2V 



/ Vj,o \ 



( Vifi \ 

\Vw-iJ 
where Uj = e 2mX i . 

When the indicial equation of zq has distinct roots Aj and Xj such that A, — Aj G 
Z, there are many possibilities for the monodromy matrix relative to the Frobenius 
basis, but in any case, the matrix still consists of blocks of entries that take the 
same form as above. 

From the above discussion we see that monodromy matrices with respect to the 
local Frobenius bases are very easy to describe. Therefore, to find monodromy 
matrices uniformly with respect to a given fixed basis, it suffices to find the matrix 
of basis change between the fixed basis and the Frobenius basis at each singular- 
ity. When the differential equation is hypergeometric, this can be done using the 
(refined) standard analytic method, in which we first express the Frobenius basis 
at z = as integrals of Barnes-Mellin type and then use contour integration to 
obtain the analytic continuation to a neighborhood of z = oo. This gives us the 
monodromy matrices around z — and z = oo. Since the monodromy group is 
generated by these two matrices, the group is determined. 

When the differential equation is not hypergeometric, we are unable to determine 
the matrices of basis change precisely. To obtain the matrices numerically we use the 
following idea. Let z\ and z-i be two singularities and {yi} and {yj}, i,j = 1, . . . , n, 
be their Frobenius bases. Observe that if yi — anyi + • • • + a in y n , then we have 

I v\ lA ■■■ v) \ /an ai9 ••• ai„\ / yi ■■■ 



2/1 

V2 



\Vn 



2/2 



(n 

2/i 

(n- 
2/2 



V'n 



(n 
Vn 



( an 

021 



a-22 



01, 
02i 



\o n i a„2 



yi 

2/2 



yi 

2/2 



2/1 

~{n- 
2/2 



1) 



\y n y' r , 

(fe) 



_(n 
Vn 



7 



Thus, to determine the matrix (a^ ) it suffices to evaluate y\ K> and y^ at a common 
point. To do it numerically, we expand the Frobenius bases into power series and 
assume that the domains of convergence for the power series have a common point 
Zq. We then truncate and evaluate the series at Zq. This gives us approximation of 
the matrices of basis changes. We will discuss some practical issues of this method 
in Section 5. 



4. The hypergeometric cases 



Throughout this section, we fix the branch cut of log z to be (— oo, 0] so that the 
argument of a complex variable z is between — ir and ir. 
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Recall that a hypergeometric function p F p -i(ai, . . . , a p ; 0i, . . . , @ p -i; z) is de- 
fined for ^ ^ 0, -1, -2, ... by 

r i : . - , - \ (ai) n ■ ■ ■ {a P )n 

p-Tp — 1 



„=S(l)nC9l)n...(/3p-l)r 



. .. ,a p ;/3i, . . . ,j3 p -v,z) = V] tt^; tttt 

(.J-JnlPlJn 

where 

(a)n 



lo(a+l)...(a + n-l), if n > 0, 
jl, ifn=0. 
It satisfies the differential equation 

(11) [0(9 + fa - 1) . . . (6 + &_ x - 1) - «(0 + a x ) . . . (0 + ap)]/ = 0. 
Moreover, it has an integral representation 

i r(A)---r(/3 p -i) y r( s + a 1 )...r(s + a p ) 
27ri r(ai)...r(ap) y c r(s + /3x)...r(s + /3 p _x) 1 n ' 

for | arg(— z)| < 7r, where C is any path from — zoo to zoo such that the poles of 
r(— s) lie on the right of C and the poles of T(s + ajt) lie on the left of C. (See [2~Tl 
Chapter 5].) Then one can obtain the analytic continuation of p F p _i by moving the 
path of integration to the far left of the complex plane and counting the residues 
arising from the process. It turns out that this method can be generalized. 

Lemma 1. Let m be the number of 1 's among (3k- Set 

F(h z) = V («i + h) n . . . (ap + h) n +h 
1 ^ ) (l + h) n (/3 1 + h) n ...((3 p - 1 + h) n 

Then, for j = 0, . . . , m, the functions 

dtp h=0 

are solutions of (!!])■ Moreover, if | arg(— z) \ < it and h is a small quantity such that 
cik + h are not zero or negative integers, then F(h, z) has the integral representation 

z h r(fli + h)... r(/3 p _x + h)T(l + h) 



' ^ 2m T(ai + h)...T(a p + h) 

f r(s + ax + h) . . . T(s + a p + h) n 

-[— z) as, 



l c T(s + 0i+h)... r(s + P P -i + h)T(s + l + h) sin7rs 

where C is any path from —zoo to zoo such that the integers 0, 1, 2, . . . lies on the 
right of C and the poles ofT(s + + h) lie on the left of C. 

Proof. The first part of the lemma is just a specialization of the Frobenius method 
(see [2]) to the hypergeometric cases. We have 



0(0 + (6+ <3p-i - l)F(h, z) = h(h + 0!-l)...(h + /Vi - 

+ (ai +h) n ... (ap + h) ^ n+h 

£j (1 + h) n -i(0i + h) n -i . . . (/3 p _i + h) n - 



-1 
and 



z{9 + ai)...(0 + a p )F(h, z) = ]T - 



(ai + h) n+ i . . . (otp + h) n+ i n+ i +h 

n~t) ( X + h )n(Pl +h) n ... (0 p -l + K) n Z 
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It follows that 

[9(9 + (3 1 -1)...{9 + f3 p -i - 1) - z{9 + ol x )...{6 + a p )]F{h, z) 

= h(h + ft - 1) . . .(h + f3 p -i - l)z h . 

If the number of l's among [3k is to, then the first non- vanishing term of the Taylor 
expansion in h of the last expression is h rn+1 . Consequently, we see that 

oh? h=0 
are solutions of (fTTj) for j = 0, . . . , to. 

The proof of the second part about the integral representation is standard. We 
refer the reader to Chapter 5 of [21]. □ 

We now prove Theorem 1. Here we will only discuss the cases 

(A,B) = (1/2,1/2), (1/3,1/3), (1/4,1/2), and (1/6,1/3), 

representing the four classes whose indicial equations at z — oo have one root with 
multiplicity 4, two distinct roots, each of which has multiplicity 2, one repeated 
root and two other distinct roots, and four distinct roots, respectively. The other 
cases can be proved in the same fashion. 

Proof of the case (A, B) = (1/6,1/3). Let h denote a small real number, and let 
F(h,z) be defined as in Lemma 1 with p = 4, a\ = 1/6, «2 = 1/3, «3 = 2/3, 
a 4 = 5/6, and [3k — 1 for all k. Then, by Lemma 1, the functions 

i gj 

Vj(z) = -p,djJ^ hF ^ Cz ))' J = 0, ... ,3, 

are solutions of 

9 4 - 11664z(0 + l/6)(0 + 1/3) (61 + 2/3) (6> + 5/6), 

where C = 11664. In fact, by considering the contribution of the first term, we see 
that these four functions make up the Frobenius basis at z = 0. 

We now express C~ h F(h,Cz) using Lemma 1. By the Gauss multiplication 
theorem we have 

r(s + i/6)r( s + i/3)r( s + 2/3)r( s + 5/6) 

IlLi r ( s + fc / 6 ) _ (27r) 5 / 2 6- 1 / 2 - 6 "r(6s + l) 
~ r(s + l/2)r(s+l) ~ (27r) 1 /22-i/2-2, r ( 2s + 1)' 

Thus, restricting z to the lower half-plane — ir < argz < 0, by Lemma 1, we may 
write 

z h F(l + /i) 4 r(l + 2h) 



C- h F{h,Cz) 



2iri T(l + 6h) 

r(6s + l + 6/i) 



T(s + 1 + h) 4 T(2s + l + 2/i) sin tts ' 

where C is the vertical line Res = —1/12. Now move the line of integration to 
Res = —13/12. This is justified by the fact that the integrand tends to as Ims 
tends to infinity. The integrand has four simple poles s = — n/6 — h, n = 1, 2, 4, 5, 
between these two lines. The residues are 



QT{n) T(l - n/6) 4 F(l - ra/3) sin7r(n/6 + h) ' 



i./6-h 
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Thus, we see that the analytic continuation of C h F(h,z) to \z\ > 1 with — n < 
argz < is given by 

C- h F{h, z)= a n B n (h)z- n/6 + (higher order terms in 1/z), 

n=l,2,4,5 

where 



a n = 



(-l) n ne- win / 6 



B n {h) = 



r(i + /i) 4 r(! + 2h )e 



— nth 



6r»r(l-n/6) 4 r(l-n/3)' y ' r(l + 6/i)sinir(n/6 + /i) ' 

On the other hand, since the local exponents at z = oo are 1/6, 1/3, 2/3, and 5/6, 
the Frobenius basis at z = oo consists of 

y n (z) = z-^ e g n (l/z), n=l,2,4,5, 

where g n = 1 + • • • are functions holomorphic at 0. It follows that for z with 

— 7T < argz < 



Set fj(z) = yj(z)/(2niY for j = 0, . . . , 3 and /„ = a n y n / sin(n7r/6) for n = 1, 2, 4, 5. 
Then using the evaluation 

r'(l) = - 7 , r"(l)= 7 2 + C(2), r'"(l) = - 7 3 -3C(2)7-2C(3), 

wc hnd 



ffs\ 






h 


= M 


h 


fl 


h 


\f») 







where 



M 



with 



/ rj-iuj/A 77 + 5\/3w 2 /36 r) + 5-\/3iw 4 / 36 



rj — iu>' 



74 \ 



-5/12 -iu/2 l/4-iw 2 /2V3 l/4-iw 4 /2\/3 -5/12 -iuj 5 /2 



ILU 
1 



^ 2 /V3 
1 



lo = e" /6 , 



^ 4 /V3 
1 



1 



/ 



68C(3) 



(2™) 3 ' 

Now let P be the path traveling along the real axis with argz = — n+ from z = —2 
to — oo and then coming back along the real axis with argz = n— to z = — 2. The 
monodromy effect on y n {z) = z~ n / 6 g n (l/ ' z) is 

Vn{z) — > y n {e 2 * l z) = e- 2 * in ' 6 y n {z). 

Therefore, the matrix representation of the monodromy along P relative to the 
ordered basis {/ 3 , f 2 , /i, /o} is 

^ 

uj- a 
lu 4 
\ w 2 / 

Now the path P is equivalent to that of circling once around z = 1/C and then 
once around z = 0, both in the counterclockwise direction. Therefore, if we denote 



T 00 =M 



M~ 
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by To and T\/c the monodromy matrices relative the basis {/3, /2, /i, /o} around 
z = and z = 1/C, respectively, then we have 

Too — Ti/qTq. 

Since To is easily seen to be 



To 



(I 


1 


1/2 


l/6\ 





1 


1 


1/2 








1 


1 










1 / 



we find 



T l/C = M 











°\ 
























V o 








w 



/l + a a6/d a 2 /d \ 



M _1 T Q -1 




\ -d 



1 -6 2 /d -ab/d 
1 
-6 1 - a/ 



where 



204 



C(3), 



b = 



7 



d = 3. 



(27ri) 3 v ~" _ 4' 
Comparing these numbers with the invariants, we find the matrix Txjc indeed 
takes the form (Q specified in the statement of Theorem 1. This proves the case 
(A,B) = (1/6,1/3). □ 

Proof of the case (A, B) = (1/4, 1/2). Apply Lemma 1 with p = 4, a.\ = 1/4, a 2 = 
3/4, a 3 = a 4 = 1/2, /3 fc = 1 for all k, and set C = 1024. Then 



yj{^) = J^{C- h F(h,Cz)), j = 0,.. 

form the Frobenius basis for 

6> 4 - 1024^(0 +l/4)(0 + 3/4) (0 + 1/2) 2 . 
Assuming that — ir < argz < 0, we have 



,3, 



C-' l F(h, Cz) = - 



T{l + hf 



2m T(l + 2h)T(l + Ah) 

r(4s + 1 + Ah)T(2s + 1 + 2h) n 



-e ms z s ds. 



T(s + l + h) 6 si~„.. 

where C is the vertical line Res = —1/8. The integrand has simple poles at — k — 
h-l/A and -k-h-3/A, and double poles at -k - h - 1/2 for k = 0, 1,2,.. .. The 
residues at s = —h — n/4, n = 1, 3, are a n C n {h)z^ h ^ n ^, where 



(_ 1 )(«+i)/2 7rr ( 1 / 2 ) e 



-7rm/4 



—■nih 



4r(l-n/4) 6 ' 
At s = —h — 1/2 we have 

T(4s + 1 + Ah)T{2s + 1 + 2h) 



C n {h) 



sin7r(/i + ro/4) 



r(s + l + /i) 6 
1 



8r(i/2) 6 

sin7rs 



, , , n 31og2 + 1 
(s + /i + l/2)~ 2 f ; Ne 

v 7 y 2r(i/2) 6 

7r , sin7r/i 



(s + /l + 1/2)" 1 + 



cos irh 



cos 2 7f/l 



(s + h+ 1/2) + 
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and 

e™z s = z -V*-h e -*i(h+i/*) (1 + (ttj + log z )( a + ft + 1/2) + • • • ). 
Thus, the residue at s = —ft — 1/2 is 



-iri(/i+l/2) 



Set 



«2 — 



7re 



8r(l/2) 6 cos7r/i 

7Ti/2 



7ri + log z + 12 log 2 + 4 - 7r 



sin7rft 

C0S7rft 



-h-l/2 



-irih 



8r(i/2) 6 ' 

We find 

C- h F(h, Cz) 



C 2 (ft) = 



COS 71"ft ' 



C* (ft) = C 2 (ft) (wi + 12 log 2 + 4 - 7T tan wh) . 



-a x B x (h) Z - x ^ - a 2 B 2 (h)z- 1/2 log z - a 2 B^{h)'^ 1/2 
— a 3 B 3 (h)z~ 3 / 4 + (higher order terms in l/z), 



where 
B n (h) 



r(i + ft) 



6 



-C n (h), 



B* 2 (h) = 



r(i + ft) e 



r(i + 2ft)r(i + 4ft) ny " zy ' r(i + 2ft)r(i + 4ft) 

Let Vj(z), j = 0, . . . , 3, be the Frobenius basis at z = 0, and 



cm- 



-3/4 



y 1 (z) = z- 1 / 4 (l + •••), 
j/ 2 *(z) = z-V^i + ...), Uz)= (logz + ff (l/z))y*(z) 

be the Frobenius basis at oo, where g(t) is a function holomorphic and vanishing 
at t = 0. Set /j(z) = y j (z)/(2ni) j for j = 0, . . . ,3, /„(z) = -a n y n (z)/ shi7r(n/4) 
for n = 1, 2, 3, and /|(z) = —022/2(2)- Using the fact that 



we find 
/7s\ 

h 

w 

where 



1 r)J 

%W = - jM j(C-' l F(ft,Cz)), 



/r/+(l-i)/48 77-5/48 -5/j/12 + 7m? + Afi-q ry+(l+i)/48\ 

(l-6z)/24 7/24 7^/24 + 7^/6 (1 + 6i)/24 

(i-l)/2 -1/2 -2/j -(i + l)/2 

\ 1 1 7T7 + 4/J 1 y 



//i\ 

/I 

v/ 3 y 



/i = 31og2 + 1, 



n 



22C(3) 



(2tt«) 3 ' 

Let P be the path from z = — 1 with argument — 7r to — oo and then back to 
z = — 1 with argument 7r. The monodromy matrix for P relative to the ordered 
basis {fij 2 ,f 2 ,h} is 

(-i 0\ 
-1 -2m 
0-10 
\ i) 
Thus, the matrix with respect to the ordered basis {/3, f 2 , /i, /o} is 

/l — 877 1-877 1/2-1977/3 1/6 - 11t7/3 + 8t7 2 \ 

-7/3 -4/3 -61/72 -41/72 + 7t?/3 

1 1 

\ -8 -8 -19/3 -8/3 + 877 J 



T — 

1 DC, 
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T n = 



Finally, it is easy to see that the monodromy around z = with respective to 

{/3,/2,/l,/o} is 

(1 1 1/2 l/6\ 

1 1 1/2 

1 1 

\0 1 / 

and after a short computation we find that the monodromy T^/c = T^T^ 1 around 

z = 1/C indeed takes the form claimed in the statement of Theorem 1. □ 

Proof of the case (A,B) = (1/3,1/3). Let z be a complex number with —ir < argz < 
0. By the same argument as before. We find that the Frobenius basis {yj} at z = 
can be expressed as 

y^) = j^{c- h F(h,Cz)), 

where C = 729 and 

v ' ' 2m T(l + 3h) 2 J c T(s + l + h) e sinns 

Here h is assumed to be a real number and C denotes the vertical line Re s — —1/6. 
Set 

a " = ~9r(l-n/3)6' n=1 » 2 - 
The residues at z = — 1/3 — h and z = —2/3 — h are 

ai {m + log z + 9 log 3 - ttV3 + cot tt(1/3 + h))z- 1/3 - h e~ nih 

and 

a 2 {iri + log z + 9 log 3 + ttV3 + 6 + cot tt(2/3 + h))z- 2/3 - h e~' Kih , 
respectively. Let 

B n (h) = v ' 2 - , j^—tt , n = 1, 2, 

T(l + 3n) z sm7r(n/3 + /i) 

and 

(/i) = B 1 (h) (iri + 9 log 3 - nV3 + tt cot tt (1/3 + h) ) , 
B 2 (h) = B 2 (h) {m + 9 log 3 + nV3 + 6 + tt cot tt(2/3 + h)) . 
Then we have 

2 

C- h F{h, Cz) = -J2 a n (B n (h)z- n / 3 logz + B* n (h)z^ n/3 ) + (higher order terms). 

n=l 

Now the Taylor expansions of B n (h) and B^(n) are 

• tt „ / /i \ , iw. / 1 \ / \ ,, 

sm 3 Bl UJ= 1 + ^H^ + W /l Hl27! + r? j /l + ---' 

. 7r f h \ iuj 2 , / iw 2 1 \ , 9 / iw 2 \ , 

sin -_B 2 tt— = 1 + — ^ - — ?= + tt: ^ + ?= + 1 'j H , 

3 2 \2mJ V3 V2V3 12/ V12V3 J 

■ k ( h \ ( 2hlo\ ( iu\tjj 2itiuj\ , / iu\u u,\ iru; \ 

s '° 5 B 'Wi) " r + IS) + (tV + — ) " - (if + W + i^J " 2 
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3 2 V 2ni 





2ttw 2 
+ V3 


+ 6^) 


V V3 


2vTlW 

3 




( ifi 2 ^ 2 


+ ^ 


7TW 2 2 

2V3 






V 2^3 


12 


-1) 


+ 


^2V + 


12V3 


27T77CJ 2 
+ 7=~ + 


niu) 2 
6 



sin (—-) = ( M2 + +6 ] + ( -f^= — + 2zw 2 V3 ] /i 

/i 2 



+ 677- 



• 9 

2V3 



where 

u = e 7Tl/3 , r, = ]!^l, m = 91og3-vrV3, M 2 = 9 log 3 + tt \/3. 
From these we can deduce the matrix of basis change between the Frobenius basis 
fj(z)= , n \. , -^C- h F(h,Cz), 7 = 0,..., 3 

and the basis 

fi(z) = -a.z-^il + •••), ft(z) = -a 2 ^ 2/3 (l + ■■■). 
A - (logz + . 9l (l/z))/*(z), / 2 - (logz + 52 (l/z))/*(z), 

where <7i(t) and c/ 2 (t) are functions holomorphic and vanishing at t = 0. The 
monodromy matrix around oo with respect to the ordered basis {/i, /* , / 2 , / 2 } is 
easily seen to be 

/ e - 27 ™/ 3 27rie- 27 ™/ 3 \ 

e -2W3 o 1 

e 27 ™/ 3 27rie 2 "/ 3 

\ e 2 "/ 3 / 

By the same argument as before, we find that the monodromy matrix with respect 
to the basis {/3, / 2 , /i, /o} indeed takes the form claimed in the statement. This 
proves the case (A, B) = (1/3, 1/3). □ 

Proof of the case (A, B) = (1/2, 1/2). Let z be a complex number such that — n < 
argz < 0. We find that the Frobenius basis {yj} at z = can be expressed as 

Vj(z) = ^(C- h F(h,Cz)), 

where C — 256 and 

-hr?,u^„\- zH r (! + ^) 8 /" T(2s + 1 + 2/i) 4 7T 



Cz ^ = -^wrwL hj+i+W ^ e " v ds - 

Here h is assumed to be a small real number and C denotes the vertical line Re s = 
— 1/4. The integrand has quadruple poles at s = —k — 1/2 — h for non-positive 
integers k. Moving the line of integration to Re s = —3/4 and computing the residue 
at s — —1/2 — h, we see that 

C- h F(h,Cz) = oi V ^^z- 1 / 2 (logz)" + (higher order terms in 1/z), 



n=0 



where 
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and 



B , ( 1 1 ) = B:;[ 1 1 ) j - ^7r 2 + ^ — ir[i tan irh + tt 2 sec 2 irh 



3 

B Q (h) = B z (h)(^- - ^/i 2 tan7r/i+ (sec 2 irh - 7/6)ttV 
+ (5/6 - sec 2 irh)ir 3 tanTr/i + 8C(3)) , 



where /i = 16 log 2 + 7ri. Let 
/o(z) = z- 1 / 2 (l + - 



hi?) 



Hz) 
i 



/i(z) = — (logz + 5l (l/z))/ (z), 
2m 

r (log 2 z/2 + 5l (l/z) logz + 52 (l/z))/ (z), 



(2tt0 2 

(log 3 z/6 + 9l (l/z) log 2 z/2 + g 2 (l/z) logz + 53 (l/^)/o(^) 



(2ttz) 3 

be the Frobenius basis at z — oo with g n (0) = 0. Using the evaluation 

r'(i) = - 7 , r"(i) = ^ + ^, r"'(i) 1 



7_ 
6 



7T 

12 ' 2 ' 3 ,v "' 12 

we can find the analytic continuation of the Frobenius at z = in terms of f n (z). 
Now the monodromy around oo relative to the basis {/^(z), / 2 (z), fi(z), fo(z)} is 

/-I -1 -1/2 -l/6\ 

-1 -1 -1/2 

0-1 -1 

\ -1 J 

From this we can determine the monodromy matrix around z = l/C with respect to 
the Frobenius basis at z = 0. We find that the result agrees with the general pattern 
depicted in Theorem 1, although the detailed computation is too complicated to be 
presented here. □ 

Of course, there is no reason why our approach should be applicable only to 
order 4 cases. Consider the hypergeometric differential equations of the form 

(12) L : 6 5 - z(9 + 1/2) (0 + A) {9 + 1-A)(6 + B){6 + 1-B). 

The cases (A, B) = (1/2,1/2), (1/4,1/2), (1/6,1/4), (1/4,1/3), (1/6,1/3), and 
(1/8,3/8) have been used by Guillera [12J, [13] to construct series representations 
for l/7r 2 . Applying the above method, we determine the monodromy of these 
differential equations in the following theorem whose proof will be omitted. 



Theorem 3. Let L be one of the differential equations in Uty) . Let yi, i = 0, . . . , 4, 
be the Frobenius basis at 0. Then the monodromy matrices around z = and z — 
l/C with respect to the ordered basis {y4/(27ri) , ?/3/(27ri) 3 , T/2/ (27rz) 2 , yi/(2ivi), yo} 
are 



(I 


1 


1/2 


1/6 


l/24\ 




f 


a 2 





—ab 


(1 - a 2 )x 


-6 2 /2 \ 





1 


1 


1/2 


1/6 






-c 2 x/2 


1 


—acx 


c 2 x 2 /2 


-(1 -a 2 )x 








1 


1 


1/2 


1 




—ac 





I -2a 2 


acx 


—ab 











1 


1 















1 





\p 











1 J 




V 


-c 2 /2 





—ac 


c 2 a;/2 


a 2 j 
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respectively, where x is an integer multiple o/£(3)/(27ri) 3 , a and c are positive real 
numbers such that a 2 , ac, and c 2 are rational numbers, and b is a real number 
satisfying a 2 + be — 1. The exact values of a, c, and x' = (2iri) 3 x/C(3) are given 
in the following table. 



A B 


a 2 


c 2 


x' 


1/2 1/2 


25/36 


64 


10 


1/2 1/4 


8/9 


32 


24 


1/4 1/6 


289/288 


8 


80 


1/3 1/4 


27/32 


24 


28 


1/3 1/6 


75/64 


12 


70 


1/8 3/8 


529/288 


8 


150 



5. Differential equations of Calabi-Yau type 

The Picard-Fuchs differential equations for families of Calabi-Yau threefolds 
known in literature have the common features that 

(a) the singular points are all regular, 

(b) the indicial equation at z — has as its only solution, 

(c) the indicial equation at one of the singularities has solutions 0, 1, 1, 2, cor- 
responding to a conifold singularity, 

(d) the unique holomorphic solution y around with y(0) = 1 has integral 
coefficients in its power series expansion, 

(e) the solutions Ai < X 2 < A3 < A4 of the indicial equation at t = 00 are 
positive rational numbers and satisfy Ai + A4 = A2 + A3 = r for some r G Q, 
and the characteristic polynomial of the monodromy around t = 00 is a 
product of cyclotomic polynomials. 

(f) the coefficients ri(z) of the differential equation satisfy 

11 3 1 

n = ^r 2 r 3 - -r 3 + r' 2 - -r' 3 r 3 - -r 3 ', 

(g) the instanton numbers are integers. 

In pQ a fourth order linear differential equation satisfying all conditions except (c) 
is said to be of Calabi-Yau type. Using various techniques, Almkvist and etc. found 
more than 300 such equations. (See Section 5 of [2] for an overview of strategies 
of finding Calabi-Yau equations. The paper also contains a "superseeker" that 
tabulates the known Calabi-Yau equations, sorted according to the instanton num- 
bers.) Among them, there are 178 equations that have singularities with exponents 
0,1,1,2. It is speculated that all such equations should have geometric origin. 

In |26j van Enckevort and van Straten numerically determined the monodromy 
for these 178 equations. They were able to find rational bases for 145 of them, 
among which there are 64 cases that are integral. Their method goes as follows. 
Let Z\ , . . . , Zk be the singularities of a Calabi-Yau differential equation. They first 
chose a reference point p and piecewise linear loops each of which starts from p and 
encircles exactly one of Zi. Then the problem of determining analytic continuation 
becomes that of solving several initial value problems in sequences. This was done 
numerically using the dsolve function in Maple. Then they used the crucial ob- 
servation that the Jordan form for the monodromy around a conifold singularity 
is unipotent of index one to find a rational basis. Finally, assuming that ([5]) and 
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dHJ hold for general differential equations of Calabi-Yau type, conjectural values of 
geometric invariants can be read off. 

Here we present a different method of computing monodromy based on the ap- 
proach described in Section 3. Let = zq, Zi, . . . , z n be the singular points of a 
Calabi-Yau differential equation, and assume that f^, i = 0, . . . , n, k = 1, . . . , 4 
form the Frobenius bases at Zi. According to Section 3, to find the matrix of basis 
change between {/j,fc} and {fj t k}, we only need to evaluate and fj™' at a 

common point £ where the power series expansions of the functions involved all 
converge. In practice, the choice of £ is important in order to achieve required 
precision in a reasonable amount of time. 

Let Ri denote the radius of convergence of the power series expansions of the 
Frobenius basis at z,. In general, Ri is equal to the distance from Zj to the nearest 
singularity zj ^ Zi, meaning that if we truncate the power series expansion of fi : k 
at the nth term, the resulting error is 



Of course, the O-constants depend on the differential equation and z,. Since we do 
not have any control over them, in practice we just choose £ in a way such that 



If this does not yield needed precision, we simply replace n by a larger integer and 
do the computation again. 

Example. Consider 



The singularities are zq = 0, Z\ = 1/3125, and zi = oo. The radii of convergence 
for the Frobenius bases at and 1/3125 are both 1/3125. Thus, to find the matrix 
of basis change, we expand the Frobenius bases, say, for 30 terms, and evaluate 
the Frobenius bases and their derivatives at £ = 1/6250. Then we use the idea 
in Section 3 to compute the monodromy matrix around z\ with respective to the 
Frobenius basis at 0. We find that the computation agrees with (J7J in Theorem 1 
up to 7 digits. 

The above method works quite well if the singularities of a differential equation 
are reasonably well spaced. However, it occurs quite often that a Calabi-Yau dif- 
ferential equation has a cluster of singular points near 0, and a couple of singular 
points that are far away. For example, consider Equation #19 

5290 4 - 23z(9210 4 + 2O460 3 + 16440 2 + 6210 + 92) 

- z 2 (38O8510 4 + 13285840 3 + 17726736* 2 + 10335280 + 221168) 

- 2z 3 (4758610 4 + 131O1720 3 + 10287916* 2 + 2089320 - 27232) 

- 68z 4 (88730 4 + 14O2O0 3 + 51390 2 - 16640 - 976) 
+ 6936z 5 (0 + 1) 2 (30 + 2)(30 + 4). 

The singularities are z = 0, z\ = 1/54, z 2 = (11 - 5\/5)/2 = -0.090..., z 3 = 
—23/34, and Z4 = (11 + 5y/E)/2 = 11.09 .... In order to determine the monodromy 
matrix around Z4, we need to compute the matrix of basis change between the 




\C- Zl \ _ \(- Zj \ 

Ri Rj 



4 - 5(50 + 1)(50 + 2)(50 + 3)(50 + 4). 
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Frobenius basis at 1/54 and that at Z4. The radius of convergence for the Frobenius 
basis at 1/54 is 1/54, while that at 2:4 is Z4 — 1/54 = 11.07 . . .. Even if we choose 
£ optimally, we still need to expand the Frobenius bases for thousands of terms 
in order to achieve a precision of a few digits. In such situations, we can choose 
several points lying between the two singularities, compute bases for each of them, 
and then use the same idea as before to determine the matrices of basis change. 

Take Equation 19 above as an example. We choose Wk = (1 + 3 fe )/54 and 
£fc = (1 + 3 /2)/54 for k = 0, ... , 5. The radius of convergence for the basis at Wk 
is 3' c /54. Thus, evaluating the first n terms of the power series expansions at C,k 
and Cfe+i will result in an error of 

Oe((l/2 + e)"), 

which is good enough in practice. 

Using the above ideas we computed the monodromy groups of the differential 
equations of Calabi-Yau type that have at least one conifold singularity^ Our result 
shows that if a differential equation comes from geometry, then the monodromy 
matrix around one of the conifold singularities with respect to the Frobenius basis 
at the origin takes the form ([7|). We then conjugate the monodromy matrices by the 
matrix (|9|) and find that the other matrices are also in Sp(4,Q). We now tabulate 
the results for the equations coming from geometry in the following table. Note 
that the notations T(d\ 1 d2) and r(di, d 2 , ^3), d2,da\di, represent the congruence 
subgroups 



T(d u d 2 )= <NeSp(4,Z) :7 = 



p| J 7 eSp(4,Z) 





* 


* 









* 


* 


# 


1 

mod di 





1 












* 


* 







and 



T(d 1 ,d 2 ,d 3 ) = {(ay) S Sp(4,( 



f 1 









V 



mod di 



a 



eZV(i,j)^(l,3), a 13 € — z, 

da 



0-21, Cl31 : 0,4,1, ^32, «34 
042 = 0, 022, 044 = 



= mod di, 
1 mod ffe, 



an, 033 



1 



A dl \ 

mod — >. 
da J 



Note also that since the matrix 



(I 





1 



V 



We have written a Maple program for the computation. It is available upon 
request. 
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is always in the monodromy groups, it is not listed in the table. The reader should 
be mindful of this omission. 



# 


H 3 


C2-H 


C3 


Generators 


in 


Rcf 










' 1 i o o N 




1 


-2 


3 


1/2 - 


A 








15 


18 


72 


-162 


10 
18 18 10 
\0 -9 -1 lj 




V 


(i 

-18 
18 


-5 
18 
-18 


-1 
4 

-3 


2 

6 

' J 




r(6,3,2) 


0! 










' 1 1 






/ -5 


2 


1/4 












16 


A O 

48 


96 


1 o o 

— 128 


10 
48 48 1 






24 -7 
-144 48 


-1 
7 




4 

24 




T(24,8,4) 


m 










\0 -16 -1 1 


{ 




V 48 


-16 


-2 


9 , 


1 














' 1 1 C 






(-11 


3 


1/3 


1 \ 
















10 




-54 


10 


1 


3 
















90 90 1 




-972 162 


19 


54 








17 


90 


108 


-90 


^0 -24 -1 l) 
'-11 3 1/3 




V 162 
-l\ 


-27 


-3 








r(18,6,3) 


01 










36 -8 -1 




3 






















-432 108 13 




-36 






















\ 108 -27 -3 




t4 






















' 1 1 


\ 






4 


1/2 


-l\ 








1 O 

18 


A A 

40 


o o 

88 


— 128 


10 
40 40 1 






12 
-72 


-7 
48 


-1 
7 


2 
-12 




r(4,2,2) 


m 

m 










\0 -14 -1 1 


\ 




V 48 


-32 


-4 


9 


A 
















1 1 1 






/-6 


4 


1/2 




















10 




14 


-7 


-1 


2 


















46 46 1 




-98 


56 


8 


-14 








19 


46 


88 


-106 


\0 -15 -1 1 
1 -45 12 J 
138 -35 - 


J 
(> 




\ 56 
-6 
18 


-32 

\ 


-4 


9 


J 




r(2,2,2) 


01 










-1058 276 47 


-13 


8 




















^ 276 -72 -12 


37 


/ 




















' 1 1 


\ 




( 7 


-1 


-1/6 












20 


54 


72 


— 18 


10 
54 54 1 






-6 
126 


1 

-18 




-2 


-2 
24 




r(6, 3, 6) 


01 










\0 -15 -1 1 


\ 




V-36 


6 


1 




















' 1 1 C 






(-11 


5 


1/2 




-i 


\ 














1 




24 


-9 


-1 




2 
















80 80 1 




-288 120 


13 




24 








21 


80 


104 


-88 


\0 -22 -1 1 
'-19 4 1/2 


J 




\ 120 
-2\ 


-50 


-5 


11 


/ 


r(8,2,2) 


01 










80 -15 -2 




8 






















-800 160 21 




-80 






















^ 160 -32 -4 


























' i i or 


\ 









\ 


















10 




10 


1 


2 


















70 70 1 




-50 


1 


-10 










22 


70 


100 


-100 


\0 -20 -1 l) 
' -9 5 1/2 




V o 

-l\ 





i ; 






r(io,io,2) 


01 










20 -9 -1 




2 






















-200 100 11 




-20 






















V 100 -50 -5 




11 / 






















' 1 1 


\ 




/ 9 


-1 


-1/8 




1 \ 








23 


96 


96 


-32 


10 
96 96 1 






-8 
288 


1 

-32 




-3 


-2 
40 




T(8,8,8) 


0! 










\0 -24 -1 1 


1 




V-64 


8 


1 




-V 
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24 



15 



66 



-150 



110 
10 
15 15 10 
-8-11 



12 -9 

-36 30 
30 -25 



-2 

4 
-12 

11 



r(3,i) 



25 



20 



us 



-120 



110 

10 

20 20 1 

-9-1 V 



-7 5 

16 -9 

-64 40 

40 -25 



-2 

4 
-Hi 

11 



r(4,i) 



3 



2G 



28 



76 



-116 



110 
10 
28 28 10 



-9 
20 



6 
-11 



-100 


(id 




11 


60 


-36 


-6 


1 











42 


1 





9 


-196 





1 


-42 











1 



-2 
4 
-20 
111 



r(4,i) 



3, 



1 


1 





1 


42 


42 



27 



42 



84 



98 



-13 


7 


1 


- 2 \ 


28 


-13 


-2 


4 


-196 


98 


15 


-28 


98 


-49 


-7 


15 ) 



r(u,7) 



3, 22 



28 



42 



84 



-96 



110 
10 
42 42 1 
V -14 -1 



-41 


12 


2 


-6 


120 


-35 


-6 


18 


-882 


252 


43 


-126 



V 252 



-72 



r(42,2) 



3, 



37 / 



186 



57 



90 



-84 



1 
10 
57 57 10 
17 -1 1/ 
-17 8 1 
36 -15 -2 
■324 144 19 
\ 144 -64 -8 



-53 
162 
1458 
\ 324 
-2 N 
4 

-36 

17 ; 



12 

-35 
324 
-72 



55 
-12 



-6 
18 
-162 
37 j 



r(3,i) 



In the second table we list a few equations whose monodromy matrices with 
respect to our bases have integers as entries. Note that the numbers H 3 , C2-H, 
and C3 are all conjectural, obtained from evaluation of the monodromy around a 
singularity of conifold type. Note, again, that the matrix 

/l 0\ 



1 
1 
\0 



1 



V 



is omitted from the table. 



29 



33 



42 



H 3 



24 



32 



c 2 -H 



72 



36 



80 



C3 



-116 



-72 



-116 



Generators 



110 
10 
24 24 1 
-10 -1 



-47 
120 
-576 
V 240 



20 4 -10 

-49 -10 25 

240 49 -120 

-100 -20 51 / 



110 
10 
6 6 10 
V> -4 -1 1, 



10 
2 10 2 
-2 1-2 
1 



1 1 

10 

32 32 1 

-12 -1 



-15 


6 


1 


48 


-17 


-3 


-256 


96 


17 


96 


-36 


-6 



-3 

9 
-48 

10 



r(24,2) 



r(2,2) 



r(ie,4) 
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51 



10 



64 



-200 



i o o 

1 

10 10 1 

\0 -7 -1 1/ 



-3 5 1-2 

8 -9-2 4 

16 20 5 -S 

V 20 -25 -5 11 



r(2,i) 



63 



62 



-310 



i o o 

10 
5 10 
-6 -1 1 



77 



-1 5 

4 -9 

-4 10 

\10 -25 



-2 
4 

-4 
11 



r(i,i) 



73 



30 



12 



1 o o 

1 

9 9 10 

V0 -4 -1 1 



2 

3 -2 
3 

-3 



T~L 1 

1 

13 13 1 

\0 -7 -1 lj 



1 
-1 
2 3 

-5 4 1 
12 -7 -2 
-36 24 
24 -16 



r(3,i) 



99 



13 



■58 



-120 



1 -2 

-2 4 

7 -12 

-4 9 . 



r(i,i) 



100 



36 



72 



-72 



1 1 oo 

1 

36 36 1 

\0 -12 -1 1/ 



1 

12 
-36 

V o 




4 
-12 
1 



24 -11 -2 4 
144 72 13 -24 
\ 72 -36 -6 13 



r(i2,i2) 



101 



25 



70 



-100 



1 1 I 

1 I 

25 25 
\0 -10 - 
1 








V 

\ 
16 



-19 

40 
200 
\ 100 



10 
-19 
100 
-50 



2 

-4 
21 
-10 



-40 
21 / 



60 1 
-225 1 -60 
1 



1 1 

1 

7 7 10 

\0 -5 -1 lj 



r(5,5) 



109 



46 



-120 



-3 


3 




1 


8 


-5 




-2 


-16 


12 




5 


12 


-9 




-3 


/ 1 










4 


1 







-4 





1 




V o 











7 J 



r(i,i) 



117 



12 



36 



-32 



1 i o 

10 
12 12 1 
\0 -5-1 
I -59 21 
-120 43 
-400 140 
\ 140 -49 







1/ 

9 18 \ 
18 36 
61 120 

-21 -W 



r(4,i) 



118 



10 



40 



-50 



110 
10 
10 10 1 
\0 - 
-19 

40 
-100 

50 



-1 

10 
-19 

50 
-25 







1/ 

4 

-8 
21 
-10 



1 
30 1 
-50 1 
V 
- 8 \ 
16 
-40 
21 





18 
-30 
1 / 



r(io,5) 



185 



36 



84 



-120 



1 


1 





o\ 


/-ll 


7 


1 





1 








24 


-13 


-2 


36 


36 


1 





-144 


84 


13 



-2 
4 

-24 

212 



r(i2,i) 
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Appendix — The index of r(di,d 2 ) in Sp(4, Z) 

by Cord Erdenberger 

Institut fur Mathematik, Universitat Hannover 

Welfengarten 1, 30167 Hannover, Germany 
E-mail : erdenberOmath . uni-hannover . de 

In this appendix we will show that the groups r(di,d 2 ) are indeed congruence 
subgroups in Sp(4, Z) and provide a formula for their index. 

Recall that for n e N the principal congruence subgroup of level n is defined by 

T(n) := { M e Sp(4, Z) | N = I 4 (mod n) } . 

It is the kernel of the map from Sp(4, Z) to Sp(4, Z/rtZ) given by reduction modulo 
n and thus a normal subgroup in Sp(4, Z). It is a well-known fact that this map is 
surjective and hence the sequence 

I 4 -» r(n) ^ Sp(4, Z) -» Sp(4, Z/nZ) -» I 4 

is exact. So the index of T(n) in Sp(4, Z) is just the order of Sp(4, Z/nZ) which is 
known to be 



[Sp(4,Z) : T(n)} = |Sp(4,Z/nZ)| = n 10 [](l ~P _ 
where the product runs over all primes p such that p|n. 
For di,d 2 e N, define 



3 )(1-P- 4 ), 



ri(di) := 



T 2 (d 2 ) :-- 



M G Sp(4, Z) : M 



M e Sp(4, Z) : M 



(I 


* 


* 










* 



* 

1 


* 



mod d\ > 


V° 


* 


* 


*) 




(i 


* 


* 




> 






1 




* 

1 


* 



mod ok / 







* 


V 





and set 



Note that 



T(d u d 2 ) :=r 1 (d 1 )nr 2 {d 2 ). 
r(di) c Fi(di) and r(d 2 ) c r 2 (d 2 ). 



Hence 

T(d) - r(di) n r(d 2 ) c f i(di) n r 2 (d 2 ) = r(d l7 d 2 ), 

where d is the least common multiple of d\ and d 2 . This shows that T(di, d 2 ) is a 
congruence subgroup, i.e. it contains a principal congruence subgroup as a normal 
subgroup of finite index. Moreover, this implies that r(di,d 2 ) has finite index in 
Sp(4, Z) and an upper bound is given by the index of T(d) as given above. 

We will from now on restrict to the case relevant to this paper, namely d 2 |di. 
Then T(di,d 2 ) is in fact a subgroup of Ti(di), namely 

/an Oi2 ai3 au\ 



r(di,d 2 ) = 



021 
031 

\a 4 i 



022 

a 32 

042 



023 

a 3 3 

043 



a 2 4 
an J 



e ri(di) 



( 022 
V«42 



024 

a 44 
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To obtain a formula for the index of this group in Sp(4, Z), we first calculate the 
index of Ti(di). Note that 

fi(di)/r(di) < s P (4,z)/r(di) ~ s P (4,z/diZ) 

and hence 

[Sp(4,Z) : ri(di)] = [Sp(4,Z/diZ) : r 1 (d 1 )/r(d 1 )]. 

The quotient ri(di)/r(di) considered as a subgroup of Sp(4, Z/diZ) via the above 
isomorphism is given by 



ri(di)/r(di) 



M e Sp(4,Z/diZ) : M 



(I 


* 


* 




> 





* 


* 


* 










1 





> 




* 


* 







An element of this group has the following form 

/l ai2 ai3 a u \ 

M a a 23 /3 
10 

\0 7 a 43 (5 / 



Let J4 := 

a 
7 <5 



® 0J' ^ C s y m Pl cct i c relation that t M J4 M = JJ 4 then implies that 
€ SL 2 (Z/diZ). Furthermore it gives rise to the following linear system: 



012 
014 



a d43 

/?043 



7023 
<5o23 



Writing this in matrix form, we have 



13 5 



a 7\ -a i3 



a 23 



012 
O14 



If we choose ai2, 013, ai4 freely, the above linear system has a unique solution 023, 043 
as ^ is in SL 2 (Z/diZ). This shows that 

|f i(di)/r(di)| = 4 ■ |SL2(z/diZ)| = 4 lld-p- 2 ) 

where the product runs over all primes p dividing di . So we have the index formula 
[Sp(4,Z) : fi(di)] = [Sp(4,Z/diZ) : ri(di)/r(di)] = d\ -p" 4 ). 

p\di 

Now we are ready to calculate the index of r(di, d 2 ) in Sp(4, Z). Since we assume 
that d2|di, we have the following chain of subgroups: 

r(di) <r(d!,d 2 ) <fi(di) <s P (4,z) 

Note that 

[fi(di) : r(di,d 2 )] = [fi(di)/r(d!) : r(rf 1 ,d 2 )/r(d 1 )] 

and by our above description this is just the index of the group 
|m g SL 2 (Z/diZ) : M = mod d 2 | 
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in SL2(Z/diZ). An easy calculation shows that this index is equal to 

p\d 2 

Putting all these together, we get 

[Sp(4,Z) : r(di,d 2 )] = [Sp(4,Z) : fi(d x )] • [Ti(di) : T(d u d 2 )] 

p\di p\d 2 

We summarize the above calculation to obtain 

Theorem. The group T(d\, d 2 ) is a congruence subgroup in Sp(4, Z) and its index 
is given by 

|Sp(4,Z) :T(d l7 d 2 )\ =df J](l-p- 4 ) ^[](l-p- 2 ). 

p|di p\d 2 

In fact, we can do a similiar calculation without the assumption that d 2 \d\ and 
obtain the same formula as given above where one has to replace d\ with the least 
common multiple of d\ and d 2 . 
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